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We argue that the ground state of a two-dimensional electron gas with Rashba spin-orbit coupling realizes one 
of several possible liquid crystalline or Wigner crystalline phases in the low-density limit, even for short-range 
repulsive electron-electron interactions (which decay with distance with a power larger than 2). Depending on 
specifics of the interactions, preferred ground-states include an anisotropic Wigner crystal with an increasingly 
anisotropic unit cell as the density decreases, a striped or electron smectic phase, and a ferromagnetic phase 
which strongly breaks the lattice point-group symmetry, i.e. exhibits nematic order. Melting of the anisotropic 
Wigner crystal or the smectic phase by thermal or quantum fluctuations gives rise to a non-magnetic nematic 
phase which preserves time-reversal symmetry. 



I. INTRODUCTION 

Enhancing the role of electron-electron interactions relative 
to that of the kinetic energy often leads to interesting many- 
body effects. Electron crystallization is an extreme example 
of this phenomenon. At low densities, where electrons are 
far apart and the kinetic energy cost of localization is low, 
Coulomb interactions dominate and the electrons form an or- 
dered state, known as a Wigner crystal 1 . 

The nature of the crystallized electronic state has been in- 
tensely investigated through Quantum Monte Carlo (QMC) 
calculations 2 3 . These calculations support the existence of the 
crystalline phase, though the density at which they find crys- 
tallization is typically significantly lower than what heuris- 
tic arguments would suggest. The Wigner crystal has also 
been sought experimentally. Despite difficulties associated 
with reaching the ultra-low density regime where crystalliza- 
tion is expected, evidence that the Wigner crystal phase may 
have been realized has been reported for experiments on ever- 
cleaner samples of two-dimensional electron gases (2DEGs) 
in semiconductor quantum wellP. 

In this paper, we study the low density limit of the 2DEG 
system with Rashba spin-orbit coupling (SOC), which is 
present whenever the 2DEG lacks inversion symmetry h This 
is the case, for example, when the 2DEG is confined in an 
asymmetric quantum well, or if it is formed at the surface of a 
three-dimensional material. As shown in Fig{T^, the resulting 
dispersion relation has an extended (highly degenerate) mini- 
mum which forms a ring in momentum space. The low -energy 
density of states exhibits a divergent van Hove singularity, 
p(e) ~ e -1 / 2 , akin to the behavior of a one-dimensional sys- 
tem (see Fig. [TJi). This is in striking contrast to the usual be- 
havior p{e) ~ const, familiar for two-dimensional systems 
without spin-orbit coupling. Thus the Rashba SOC greatly 
enhances the role of interactions relative to that of kinetic en- 
ergy in the low-density limit. 

For a system with Coulomb interactions, V(r) ~ 1/r, the 
ground state at low densities is a Wigner crystal, just as for a 
2DEG without spin orbit coupling. Remarkably, however, we 
find that with Rashba SOC, broken symmetry states appear to 
be favored over the uniform Fermi liquid (UFL) state even for 



short-range interactions, V(r) ~ l/r a with a > 2. (In partic- 
ular, note that the Coulomb interaction screened by a metallic 
gate is described by a — 3.) The instability in this case occurs 
at an electron density n for which the Fermi energy is smaller 
than an energy scale set by the SOC. 

We have investigated candidate ordered states by construct- 
ing variational wave-functions, determining the patterns of 
broken symmetry which minimize their variational energies, 
and then comparing the energy to that of the UFL. We con- 




Figure 1 : (a) Dispersion of a particle with Rashba SOC. The min- 
imum of the dispersion occurs on a ring in k space, marked by a 
red circle. The red arrows show the spin polarization of the different 
Bloch states, (b) Density of states as a function of energy, corre- 
sponding to the dispersion shown in (a). Near the band bottom, the 
density of states diverges as p(E) ~ (E)^ 1 ^ 2 . (c) Schematic phase 
diagram in the low-density limit with repulsive electron-electron in- 
teractions which decay at long distances as V ~ r~ a . For a < 2, 
the ground state is an isotropic Wigner crystal (i.e. it preserves a dis- 
crete rotational symmetry, C„ with n > 2). For a > 2, states with 
a further broken rotational symmetry are favored, (i) represents an 
anisotropic Wigner crystal with a unit cell which becomes paramet- 
rically anisotropic in the low-density limit, (ii) and (iii) represent 
snapshots of a smectic state and a ferromagnetic nematic liquid state, 
respectively. 



2 



sider the following types of broken symmetry states: 1) 
Wigner crystalline (WC) states, i.e. insulating states with only 
discrete translational symmetry corresponding to one elec- 
tron per unit cell, allowing for various possible crystal struc- 
tures corresponding to different patterns of rotation symmetry 
breaking; 2) an electron smectic state which breaks transla- 
tional symmetry in only one direction, and can be viewed as 
a partially melted version of an anisotropic WC; 3) a ferro- 
magnetic nematic state which preserves translation symmetry, 
but breaks time reversal symmetry and rotational symmetry 
- this state is invariant under time reversal followed by a ro- 
tation by 7r around the symmetry axis. Note that we refer 
to a WC as anisotropic when only a discrete 2-fold rotation 
symmetry (C2) remains unbroken. In the limit of low density, 
each of these ordered states has parametrically lower energy 
(in powers of the density n) than the UFL. This strongly sug- 
gests that the UFL is unstable at low densities. In contrast, the 
energy balance between different broken symmetry phases is 
more delicate, and may well depend on long-distance fluctu- 
ational effects that are not well captured by variational wave- 
functions; we will return to this point in the final section of 
the paper. 

The nature of the low-density instability of the UFL, and 
the origin of the strong tendency of the system to a nematic 
pattern of rotation symmetry breaking (whether or not it is ac- 
companied by other patterns of symmetry breaking) can be 
most easily seen by studying candidate WC wavefunctions. 
A schematic version of the resulting WC phase diagram is 
shown in Fig. [TJ; as a function of the exponent a. For a > 2 
(short-range interactions), the unit cell of the Wigner crystal 
becomes increasingly anisotropic, with an aspect ratio that di- 
verges in the low-density limit. This unusual behavior can 
be traced back to the form of the single-particle dispersion, 
Fig[lJ, which is strongly anisotropic in the directions perpen- 
dicular and tangential to the ring-like minimum. 

A complimentary view can be obtained by considering a 
ferromagnetic state. Because of the Rashba SOC, the ori- 
entation of the magnetization vector (which is always in- 
plane) necessarily defines a preferred nematic axis - any fer- 
romagnetic state must necessarily be nematic, although a non- 
magnetic nematic phase is possible. For example, a magnetic 
moment in the y direction implies a special role for the point 
ko in Fig. la, which defines the point on the ring of minimum 
dispersion with the largest possible value of k x . At low den- 
sity, the resulting Fermi surface forms an ellipse encircling 
this special point. As in the usual Stoner theory of ferro- 
magnetism, spin polarization lowers the interaction energy via 
the Pauli-exclusion principle, which helps electrons to avoid 
each other at short distances. However, the cost in kinetic 
energy is parametrically smaller than in a conventional FL, 
owing to the divergent density of states. The variational en- 
ergy we find for the ferromagnetic nematic state differs from 
that of the anisotropic crystal only by a numerical constant 
for 2 < a < 3, so it is not possible, on the basis of the 
present considerations, to confidently determine which (if ei- 
ther) is the preferred ground state. For a > 3, the ferromag- 
netic nematic state has parametrically lower energy than the 
anisotropic WC, suggesting that it is a better candidate ground 



Table I: Scaling of the ground state energy per particle as a func- 
tion of the electron density n in the low density limit, for the various 
candidate ground states considered in this paper: the nematic ferro- 
magnetic (FM) state, the anisotropic Wigner crystal (AWC), and the 
smectic. For a < 2, the isotropic Wigner crystal always has the 
lowest energy. 
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state. The scaling of the ground state energy as a function of 
the electron density for each type of state considered in this 
paper is listed in Table [I] 

This paper is organized as follows. The model is described 
in Sec. [IT] In Sec. Ill we explain the basic physics leading 
to anisotropic Wigner crystal formation. We analyze three 
cases: contact interactions, extended short-range interactions 
(which fall off with distance with a power which is larger 
than 2), and long-range interactions. Considerations related 
to the magnetic structure of the Wigner crystal phase are dis- 
cussed in Sec.|lV] In Sec. [V] we discuss melting the Wigner 
crystal partially to obtain a smectic state. In Sec.fVIlwe dis- 



cuss the ferromagnetic nematic state. Sec. VII presents a pro- 
posed schematic phase phase diagram for a 2DEG with SOC 
and screened Coulomb interactions. In Sec. I Villi we discuss 
possible realizations in electronic and atomic systems. Ap- 
pendix|A]presents the solution of a single Rashba particle in a 
box problem, which is crucial for the arguments regarding the 
anisotropic Wigner crystal and the smectic phases, and Ap- 
pendix [B] presents the details of the Hartree-Fock analysis of 
the ferromagnetic nematic state. 



II. MODEL 

We consider a 2DEG with Rashba SOC and repulsive 
electron-electron interactions, described by the Hamiltonian 
(in units with h = 1) 



H = 



El- 



2k 



(Vj X z) ■ (7j 



Er 



(i) 



Here, m is the electronic effective mass, ko is a parameter that 
characterizes the strength of the SOC, Eq = fcg/(2m), Bj is 
the vector of Pauli matrices which act on the spin of electron 
j, and is the (repulsive) electron-electron interaction 

potential. The Hamiltonian (JTJ is invariant under translations, 
under rotations around the z axis, and under mirror reflections 
about the x and y axes, M x and M y , but not under inversion 
r — > — r. 

Below we consider cases in which, at large inter-particle 
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separation r, the interaction potential decays as a power law: 



V(\r\) 



Vo 

j. a ' 



(2) 



We distinguish between long-range and short-range interac- 
tions, which are characterized by a < 2 and a > 2, respec- 
tively. The bare Coulomb interaction is described by a = 1, 
while screening can lead to a > 1. In particular, screening 
due to a nearby metallic gate top gate leads to a = 3. 

In the absence of interactions, V^(|r|) = 0, Eq. ([T|l yields 
the single-particle dispersion law (see Fig.[lJ): 



E(k) 



2m 



(k 2 ±2k Q k) + E , 



(3) 



where k is the electron momentum (h = 1), and k = \k\. The 
minimum kinetic energy occurs for any value of momentum 
falling on a ring of radius k = ko, with a minimal value of 
E = 0. The corresponding density of states, shown in Fig. [TJd, 
is given by 



D(E) = 




(E < E ), 
(E > E ). 



(4) 



For E > E , the density of states is independent of energy, 
just as for a usual 2DEG without spin-orbit coupling. How- 
ever, for E — > 0, the density of states diverges as \j\fE. 
Because this divergence will play a crucial role in the anal- 
ysis below, we comment briefly on its origin. The divergence 
comes from the fact that the minimum of kinetic energy is in- 
finitely degenerate, occurring everywhere on a ring in momen- 
tum space, rather than at a single point or finite set of points. 
In the presence of crystalline anisotropy (which manifests it- 
self through corrections to the effective mass approximation in 
real materials), the divergence is cut off near the band bottom 
and the density of states goes to a constant. However, as long 
as koa/ir <C 1, where a is the lattice constant (i.e. as long as 
the spin-orbit coupling is weak), the crystal field anisotropy 
terms are small and Eq. Q provides a good approximation 
down to energies of order \E \(k a) 2 above the band bottom. 



III. WIGNER CRYSTAL 

A. Instability of the Fermi liquid state 

We begin by considering the stability of the uniform (Fermi 
liquid) state. According to Eq. at low densities, the Fermi 
energy ep is 



Am 2 \E \' 



(5) 



where n is the density of electrons per unit area. Thus, the 
kinetic energy per particle in the homogeneous state is 



1 f eF 
£km = - / eD (e) de 
n Jo 



TT 2 n 2 



12m 2 \E \' 



(6) 



The potential energy per particle, on the other hand, is 

£ P ot = / d ?d7V (f- ?) (n(r>(r<)), (7) 

where fl is the total area of the system, n(r) is the local den- 
sity at position r, and (■ • • ) denotes averaging in the uniform 
(Fermi gas) state. In the low-density limit, for short-range 
interactions (a > 2 in Eq. El, e pot oc n. For long-range 
interactions, the right hand side of Eq. (|7]l diverges. Here a 
neutralizing background must be taken into account, leading 



to £ pot oc n 



a l 2 . We see that in all cases, in the low-density 
pot 3> £kin> suggesting that the uniform state is un- 



limit, e 

stable to forming some sort of order. One possibility is that 
at asymptotically low densities, the ground state is a Wigner 
crystal, as in a 2DEG with no SOC. Note, however, that here, 
in the presence of Rashba SOC, this instability occurs even for 
short-range interactions. 



B. Contact interactions 

For simplicity, we begin by considering the case of the 
shortest range interactions: repulsive "contact" interactions. 
We start with a variational wavefunction which minimizes the 
interaction energy, taking each electron to be confined to a 
rectangular box of dimensions L x x L y , with different boxes 
non-overlapping. This is a zero energy eigenstate of the po- 
tential energy operator. In order for the boxes to tile the plane, 
L x and L y are constrained by the condition 



L x Ly — 



1 



(8) 



Thus, we have a single variational parameter, the aspect ra- 
tio rj = L x /L y , which we use to minimize the kinetic energy 
of the trial state. The kinetic energy per particle in the varia- 
tional state is given by the ground state energy of a single par- 
ticle in a box with Rashba SOC. This problem is investigated, 
both analytically and numerically, in Appendix A. Surpris- 
ingly, unlike the case with no SOC, the ground state energy 
in the low-density limit is minimal for rj ^ 1. In the 77 ^> 1 
limit, we find the following expression for the ground state 
energy as a function of r\ and 



e(n,rj) 



2m 



Arf 



Bn 



(9) 



where A and B are numbers of order unity, see Eq. (A14i in 
Appendix A. Minimizing Eq. (|9]l with respect to r\, we find 
that the optimal aspect ratio rf scales as 



»?*~K*o) 5 > 

and the ground state energy per particle scales as 

e* {^)~\E \{n/kiy . 



(10) 



(11) 



Therefore, in the low-density limit, we get that the energy per 
particle of this anisotropic Wigner crystal state is parametri- 
cally lower than that of the uniform state, which scales as n. 
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Note that, consistent with our assumptions, the optimal aspect 
ratio of the unit cell in the Wigner crystal becomes parametri- 
cally large at low densities. 

The fact that the kinetic energy is minimal for an 
anisotropic box can be understood as follows. Suppose that 
the ground state wavefunction is a superposition of plane 
waves with wavevectors close to some wavevector k* of 
length ko, for which the dispersion ^ is minimal. Near k*, 
the dispersion is quadratic in the radial direction, while it is 
anomalously fiat (quartic) in the transverse direction. There- 
fore confinement in the direction perpendicular to k* is less 
costly than confinement parallel to k*, and the optimal aspect 
ratio is such that the box is long in the direction of k*, and 
short in the transverse direction. 



C. Extended short-range interactions 

We now turn to the case of extended, short-range interac- 
tions, which corresponds to 2 < a < oo. We show that in this 



case, as in the case of contact interactions, the Wigner crystal 
state is extremely anisotropic in the low-density limit. 

In the case of extended interactions, the potential energy in 
the Wigner crystal phase cannot be neglected. To estimate the 
potential energy, we consider the same variational wave func- 
tion as before, in which the particles are localized in an array 
of non-overlapping L x x L y boxes. To estimate the potential 
energy, we will replace the wavefunction of each particle by a 
constant, such that the density is uniform, n = l/(L x L y ); the 
parametric dependence of the energy on n and 77 should not 
depend on this assumption. Let us focus on the anisotropic 
limit, L x ^> L y , assuming that this is the optimal configura- 
tion. The interaction energy of a given particle with all the 
other particles is 



s v {n,r,) « 2(Ux + U 2 ), (12) 
where U\ and U2 are given by 



Ui = 



OO roc 



dx dy dx' j dy' 



00 J L 



Vn 



= Von^T]^ 1 Ci 



(13) 



and 



L x Ly Jo 



dx dy dx' / dy' 



[(a; — x') 2 + (y — y') 



[(x - x'Y + (y - y')T /2 

/^V^n^lCa + Otf-")}. (14) 



V 



C\ and C2 are dimensionless constants which depend on a. In 
the 77 3> 1 limit, we get that U\ 3> U2, and therefore we ne- 
glect the latter. Now, if we assume that 77 ~ wT s , as Eq. ( 10 1 
suggests in the case of contact interactions, we find 



(15) 



We see that, for a > 3, e v becomes negligible compared to 



in the n 







the kinetic energy £*(rf ) ~ n 4 / 3 , Eq. (11 
limit. Therefore, Eqs. ( [T0| and ( [TT) are not modified in this 
case. For a < 3, e v dominates in the low-density limit, and 
we should consider both the kinetic and potential energies, 
Eqs. ([TTJ.Q: 



Etot = £v (n, 77) + e (n, 77) 
n / . 1 Bn 



2m V ^ k. 



2 77- I + dVon^r] * . (16) 



Minimizing with respect to 7/ and keeping only the most sin- 
gular term asn->0 gives 



1 



(2m Vb 



,2/a 



n*- 1 (2<a< 3). 



(17) 



We see that, for 2 < a < 3, 77 still becomes parametrically 
large in the n — > limit. Inserting Eq. ([TTJi back into Eq. ( [To} , 
we get 



-(2mV ) 2/a n 2 ( 1 ~i) (2 



< a < 3). 



(18) 



Thus for 2 < a < 3 the anisotropic Wigner crystal has para- 
metrically lower energy per particle than the uniform state. 



D. Long-range interactions 

For a < 2 (long range interactions), the Wigner crystal in 
the low density limit has the same hexagonal (Cg) symmet- 
ric triangular structure as the classical crystalline phase which 
minimizes the potential energy. We will refer to the hexagonal 
crystal as "isotropic", as opposed to the "anisotropic" crystal 
described previously, which has a lower symmetry. To show 
that the Wigner crystal is isotropic for a < 2, we note that the 
potential energy in a classical crystal scales as 



£wc 



(19) 



If, in the Wigner crystal phase, each electron is confined to a 
region whose dimension is some fraction of the mean inter- 
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electron distance, the kinetic energy cost of forming the crys- 
tal scales as the density n, as in the case without SOCP1 (Here, 
unlike before, we assume that the region to which the electron 
is confined has an aspect ratio of order unity.) Therefore, in 
the low-density limit, crystallization yields a potential energy 
gain which overwhelms the kinetic energy cost. Thus, to first 
approximation, we may ignore the kinetic energy. The ground 
state is therefore a hexagonal Wigner crystal, and is not qual- 
itatively affected by the Rashba SOC. 



IV. MAGNETIC PROPERTIES OF THE WIGNER 
CRYSTAL 

So far, we have ignored the magnetic degrees of freedom. 
The ground state of a Rashba particle in a box is two-fold de- 
generate, according to Kramers' theorem. Correspondingly, 
the variational wavefunction we considered (in which elec- 
trons occupy non-overlapping boxes) is 2 Ar -fold degenerate, 
where N is the number of electrons. This degeneracy is lifted 
by exchange interactions. 

First, we elucidate the nature of the Kramers pair of ground 
states of a single electron in a box. Because of the spin-orbit 
coupling, these states are not eigenstates of the spin opera- 
tor a. However, in the anisotropic (77 = L x /L y 1) limit, 
there are particular linear combinations of the two ground 
states which are approximately polarized in the ±y directions, 
where y is the narrow direction of the unit cell. The expecta- 
tion values of all other spin components are small; this is true 
for any choice of basis in the ground state Hilbert space. 

To demonstrate this, we calculate the quantities 



'1 



X l(aN/3>| 2 , 



2 ^ 

a,/3=l,2 



(20) 



where \a = 1,2) are the two ground states obtained from 
the numerical solution of the particle in a box problem (see 
Appendix [A|, and s, = <j,/2 where (Xi— Xl y tZ are Pauli (spin) 
matrices. As defined, Si is the maximum expectation value of 
the spin component % in the ground state manifold spanned by 
the Kramers pair \a = 1, 2). The values of Si- Xt y >z as func- 
tions of the aspect ratio r\ are shown in Fig. [2] As 77 increases, 
S y becomes close to 1/2, while S x z — > 0. This can be un- 
derstood as a consequence of the fact that when L x ^> L y , 
the ground states contain mostly components with momenta 
close to k — ±kox, with spin polarizations close to the ±y 
directions. 

The magnetic degrees of freedom in the anisotropic Wigner 
crystal can therefore be thought of as Ising-like spins, which 
are polarized in the ±y directions. These spins are coupled by 
exchange processes, which generate N-body interactions®^. 
In addition, because of the spin-orbit coupling, Van der Waals- 
like interactions generate spin-spin terms 11 , which are not ex- 
ponentially suppressed in the Wigner crystal phase. A detailed 
estimate of these interactions is complicated, and we defer 
their analysis for later work. 




Figure 2: The quantities S x , y , z , denned in Eq. pOps. the aspect ratio 
T) = L x /Ly of a box with fixed area Q, = 75/ k . When L x /L y — 1, 
Sx = Sy', as L x I Ly increases, S y becomes close to 1 and S x — > 0. 



V. SMECTIC STATE 

The anisotropic Wigner crystal variational wavefunction 
described above assumes that translational symmetry is bro- 
ken. However, long-range quantum fluctuations can restore 
translational symmetry, either partially or fully, resulting in 
either a smectic state (which breaks translational symmetry in 
only one direction), a nematic state which breaks rotational 
symmetry but is translationally invariant, or an isotropic liq- 
uid. Since the ground state energy is mostly sensitive to short- 
range correlations, the crystal and the various liquid states can 
have close energies. Below, we demonstrate that one can write 
an explicit wavefunction describing a smectic state with the 
same parametric dependence of the ground state energy on 
density as that of the anisotropic WC. In the next section, a 
ferromagnetic nematic variational wavefunction will be de- 
scribed. 

Let us consider, for simplicity, the case of short-range (con- 
tact) interactions. To construct a wavefunction for the smec- 
tic, we consider a trial Hamiltonian in which the electrons are 
confined to move along an array of strips of width L y with 
infinitely hard walls separating different strips. The problem 
is then reduced to finding the electronic ground state of a strip 
with a linear density of nL y . The single-particle dispersion 
in the strip is derived in Appendix [A| (Eq. A12i. The dis- 
persion of the each transverse sub-band has two degenerate 
"valleys" at ± (k + 8k*), where dk* ~ 1/ (k L 2 y ). We will 
assume that only the lowest sub-band of the strip is occupied. 



(For a fixed L y , this is valid for a sufficiently small density 12 .) 
Moreover, we may ignore the "valley" degeneracy since every 
electron can be assumed to be in one of the two valleys, and 
exchange processes between the valleys are suppressed. We 
thus obtain an effective one-dimensional Hamiltonian for the 
motion along a strip: 



H = Y,( A -^) + 2^ v{x *- x > h (21) 



1,3 

where, from Eq.(A12i, A = A\/9,mk^Ly and m* oc m (Ai 



6 



is a dimensionless constant). In the n limit, the interac- 
tion becomes strong compared with the Fermi energy, and the 
electrons behave as effectively hard core particles (indepen- 
dently of their valley index). The system can be mapped onto 
a non-interacting spinless fermion problem. The ground state 
energy per particle is thus 



(T \ - Al 

S{Lv) ~ 8mk 2 L* 



(™ L v) 
6m* 



(22) 



Minimizing this expression with respect to L y , we obtain the 
ground state energy per particle of the smectic state: 



£SM 



1.2 



(23) 



The scaling of the energy of the smectic state with n is thus 
the same as that of the anisotropic Wigner crystal, Eq.dll 12 



The numerical prefactor, which cannot be determined reliably 
from such simple considerations, is therefore important in de- 
termining which of these two states is favored in the n — > 
limit. 

In the case of extended interactions which decay with an 
exponent a, one can use the same variational wavefunction 
for the smectic, in which the expectation value of the potential 
energy is finite, and then minimize the total energy over L y . 
The calculation proceeds in essentially the same way as in 
Sec. |III C| and we will not repeat the details here. The result 
is that the parametric dependence of the smectic variational 



energy on n is the same as that of the Wigner crystal, Eq.( 18 1 
for any a. 



VI. FERROMAGNETIC NEMATIC STATE 

Finally, we consider a complete melting of the anisotropic 
Wigner crystal phase, preserving its preferred orientation. 
This results in a nematic state. Similarly to the situation in 
the Wigner crystal and the smectic states, we expect that in 
the low-density limit, only states in the vicinity of two op- 
posite points ±k on the ring of minimal dispersion will be 
occupied. For simplicity, we will assume the occupation is 
limited to the vicinity of only one point on the ring, k$, which 
makes the nematic state also ferromagnetic (with an in-plane 
magnetization perpendicular to k ). Such a state is particu- 
larly easy to describe within a Hartree-Fock approximation. 
We emphasize, however, that a paramagnetic nematic state is 
also possible, although it is not easily captured by a simple 
variational wavefunction. 

The Hartree-Fock analysis of the ferromagnetic nematic 
state is straightforward, and is described in Appendix |B| At 
a sufficiently low density, a spontaneous in-plane magnetiza- 
tion develops, and the Fermi surface becomes asymmetric. At 
asymptotically low densities, the Fermi surface becomes an 
ellipse centered around one of the points on the minimal dis- 
persion ring in momentum space. The total variational ground 
state energy per particle in this limit scales with density as (see 



£fm 




(24) 



Eq. BIO i 



Comparing this to Eqs.( 11 1,( 18 i, and (23 I, we see that the 
ground state energy of the ferromagnetic state is parametri- 
cally smaller than that of either the anisotropic WC or the 
smectic states for a > 3, making it the best candidate for 
the ground state. For 2 < a < 3, the scaling of the ground 
state energy with density of all three states has the same expo- 
nent. The energies differ only by the prefactor, which cannot 
be estimated reliably within the simple variational approach 
used here. 

The ferromagnetic nematic state spontaneously breaks time 
reversal (T) and rotational symmetry (IZe) about the z axis, 
as well as the mirror symmetry, A4 y , for the plane parallel 
to the ferromagnetic moment. However, it preserves the prod- 
uct, TTZ-n, of time-reversal and rotation by ir (hence the name, 
"nematic") and reflection through the plane perpendicular to 
the moment, M. x . The latter symmetry insures that there is no 
out-of-plane component of the magnetization, and no anoma- 
lous Hall effect. Note that, even though this state carries a 
finite crystal momentum, it does not carry a finite current den- 
sity, as required by a theorem by F. Blocr) 13 " 15 ! There is, how- 
ever, a large anisotropy in the in-plane Drude weight (effective 
mass). From our Hartree-Fock state (see Appendix[B|, we find 

^ 4 " i 

that the anisotropy scales as n <* for a < 4, and as ri for 
a > 4, in the low-density limit. 

The physics behind the formation of the in-plane ferromag- 
netic state is similar to the usual Stoner picture for ferromag- 
netism: the system gains exchange energy by polarizing, at 
the expense of kinetic energy. In a system with Rashba SOC at 
low density, the exchange energy gain exceeds the kinetic en- 
ergy cost due to the high density of states. In the low-density 
limit, the Fermi surface becomes parametrically anisotropic. 
Qualitatively, the short-range correlations in this state are sim- 
ilar to those of the anisotropic Wigner crystal. This explains 
why these states are close in energy, at least for a < 3. The 
long-range correlations, however, are very different: the fer- 
romagnetic state is a fluid, whereas the Wigner crystal is an 
insulator. 



VII. PHASE DIAGRAM 

So far, we have argued that for sufficiently low density 
and for short-ranged interactions, the system breaks rotational 
invariance, going into either an anisotropic Wigner crystal, 
smectic, or a nematic state. We now discuss the global fea- 
tures of the phase diagram as a function of density and the in- 
teraction range. For concreteness, let us discuss a 2DEG with 
Rashba SOC with screened Coulomb interactions, where the 
screening is from a nearby metallic gate at a distance £ away. 
The effective electron-electron interaction is V (r) s» e 2 /nr 
for r -C £, where k is the dielectric constant of the surround- 
ing material, and V (r) ~ e 2 £ 2 /nr 3 for r>(. 

The broken symmetry state forms at a density n* at which 
various scales become comparable to each other. We define a 
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\\ l/r s ~(§/« o r 




Highly 
anisotropic phases 



Figure 3: Sketch of the phase diagram of a 2DEG with Rashba SOC, 
as a function of l/r s = ao (7m) 1 ' 2 and the screening length £, 
where ao = n/me 2 is the Bohr radius. The three regions cor- 
respond to the uniform Fermi Liquid (UFL), the isotropic Wigner 
Crystal (WC), and a phase featuring a high degree of anisotropy, 
which can be either an anisotropic Wigner crystal, a smectic, or a 
nematic. The red dashed lines correspond to r s ~ f/ao, where 
the crossover between effectively short-range (screened) and long- 
range Coulomb interactions occurs; l/r s ~ fcoao, where the SOC 
length scale becomes comparable to the inter-particle distance; and 
l/r s ~ (^/ao)" 1 / 2 , where the energies of the uniform Fermi liquid 
and the anisotropic phases become comparable, see Eq.l|25[l. 



density n* at which the energy of the broken symmetry state 
per particle is comparable to that of the uniform Fermi liquid 
state (which is dominated by Coulomb energy): 



1 



2 ^ 



K) ; 



(25) 



where e is the electron charge, and n is the dielectric con- 
stant of the host material. On the left hand side we have 
used that, for short range interactions, the potential energy 
of the uniform state scales linearly with density (assuming 
that (n*)~z ^> £), and on the right hand side we have used 



Eq.( 18 » for the energy of the broken symmetry state with 
Vo ~ e 2 £ 2 1 k. Note that the energies of all the different candi- 
date states are the same up to a numerical prefactor in the case 
a = 3; compare Eqs.(18i and (BlOi. This gives n* ~ (arj£)> 
where ao = n/me 2 is the effective Bohr radius. 

In addition, we define a density n\ at which the inter- 
electron distance is comparable to the screening length £, and 
a density n\ at which the Fermi wavevector is comparable to 
fcrj. These characteristic densities are given by n* < n\ = 
l/£ 2 and n* < = kg. The strongly anisotropic states are 
favored at densities which satisfy n < n* = min[n*, ri^, n§\. 

We can now speculate about the structure of the zero 
temperature phase diagram of a 2DEG with Rashba SOC, 
sketched in Fig. [3] as a function of the dimensionless inter- 



electron spacing r s 



and the screening length 



£. Let us consider large £, for which the Coulomb interac- 
tions are effectively unscreened. Then, a Wigner crystal with 



hexagonal symmetry forms when r s > r s c « 35. - Imagine 
that we start deep in the Wigner crystal phase, with arbitrarily 
large r s and £. Upon decreasing £, while keeping r s fixed, 
eventually we reach £ < r s ao, where the interactions are ef- 
fectively short-ranged and the kinetic energy becomes impor- 
tant. At some point along this path, we expect a phase tran- 
sition from the hexagonal Wigner crystal to one of the phases 
of lower rotational symmetry: either an anisotropic Wigner 
crystal, a smectic, or a nematic, which can also be ferromag- 
netic. Which of these phases is realized cannot be determined 
reliably on the basis of the present analysis. 

At higher densities, such that r s c < r s < l/(fc ao), the 
SOC can essentially be ignored. Then, upon decreasing £ 
from the Wigner crystal, we expect a transition to a Fermi liq- 
uid. The reentrant tip of the Wigner crystal phase originates 
from the same physical reasoning as that described in Ref[T6l 

As drawn, a sliver of UFL is shown between the isotropic 
WC and broken rotational symmetry phases. Such a region 
may or may not exist, depending on details of the numerical 
factors which are beyond the scope of the calculation here. 



VIII. POSSIBLE REALIZATIONS 

The physics described here could be relevant to electrons in 
two-dimensional heterostructures which lacks inversion sym- 
metry, such as GaAlAs quantum wells. The magnitude of 
the Rashba spin-orbit coupling in these systems, however, 
is rather small. As discussed above, a necessary condition 
for realizing the phases with broken rotational symmetry is 
r s ^ 1/koao, or equivalently, kp < fcrj; in typical GaAlAs 
quantum wells, ko/kp ~ 0.5 or less, and the characteris- 
tic energy scale of t he SO C, E is at most of the order of 
a few degrees KelvirP^EH. More promising systems are sur- 
face states of heavy metal surface alloys. For instance, the 
boundary between Bi x Pbi_ x and Ag(lll) supports a surface 
state with v ery s trong Rashba SOC, with k rs 2nm~ 1 and 
E « 0.1eVPI22lMoreover, it was demonstrated that by vary- 
ing x, the Fermi level of the surface state can be tuned to 
be lower than Eq. 21 The Coulomb interactions on the surface 
are naturally screened by the metallic bulkH Detecting the 
broken rotational symmetry on the surface poses a challenge, 
because transport measurements would be dominated by the 
bulk. One possibility is to look for signatures of anisotropy in 
the finite-frequency response, e.g. in the optical conductivity, 
assuming that the anisotropic domains can be aligned (e.g., 
by application of an in-plane magnetic field). Scanning tun- 
neling microscopy can be done on metallic surface alloyP, 
and used to detect anisotropy in the electronic structure. Fi- 
nally, magnetic spectroscopy can be used to detect the fer- 
romagnetic state, which has a large in-plane moment. Such 
measurements have recently been done^l on the conducting 
interface between LaA103 and SrTiOa, and indeed, large in- 
plane moments were found. Whether these are related to the 
mechanism described in this paper remains to be seen. 

It has been proposed 25 that similar physics can arise lightly 
doped bilayer graphene with a perpendicular electric field, in 
which the single particle dispersion has a minimum on a ring 
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in k-space, even without SOC. In bilayer graphene the single 
particle dispersion is valley and spin degenerate. The ground 
state is likely to have additional broken symmetries, lifting 
these degeneracies. 

It is interesting to note that the considerations which lead to 
broken rotational symmetry at low densities are independent 
of the particle statistics, and are thus valid for two-component 
bosons with effective (isotropic) Rashba-like spin orbit inter- 
actions. Recently, various techniques were proposed to real- 
ize effective SOC in systems of trapped ultracold atomd^U. 
The properties of such systems have been the subject of in- 
tense study ^En. Highly anisotropic phases may be accessible 
in such systems at sufficiently low densities. Indeed, it was 
found 31 that the ground state breaks rotational symmetry at 
low densities, and that the ground state energy per particle 
scales as n 4 / 3 in the limit n —> 0, consistently with our re- 
sults for the anisotropic WC and smectic phases with contact 
interactions. 



IX. CONCLUSIONS 

In the presence of strong Rashba SOC, even short-range 
electron-electron interactions become important. As a result, 
the system is expected to form a broken symmetry state at 
low enough densities. In this work, we have shown that for 
sufficiently short-range interactions, states which break rota- 
tional symmetry are favored in the low-density limit. This is a 
consequence of the fact that the single-particle dispersion has 
a minimum on a ring of finite radius in fc-space, rather than 
at a single point. This physics is not limited to the case of 
Rashba SOC; for instance, a similar situation arises in spin- 
imbalanced fermionic superfluids with no SOC, in which the 
majority-spin quasiparticles have a dispersion which is mini- 
mal near Jzf, or in bilayer graphene with a transverse electric 
fielcP. 

We believe that the variational wavefunctions and phys- 
ical arguments proposed above capture the correct scaling 
of the ground state energy, which is found to be paramet- 
rically lower than other states (e.g. a uniform Fermi liquid 
or an isotropic Wigner crystal). However, this approach is 
too crude to answer some important, more detailed questions, 
such as discriminating between the different broken symme- 
try states considered here. For sufficiently short-range inter- 
actions (which fall off with an exponent larger than 3) a ne- 
matic ferromagnetic state has a parametrically lower energy 
than all the other states considered here, and is therefore the 
best candidate for the ground state. It is not clear, at this point, 
whether a non-magnetic nematic state is a competitor or not; 
such a state is harder to capture within a simple variational 
approach. More detailed calculations will be needed to de- 
termine the phase diagram for interactions which fall off with 
distance with a power of 3 or less. For instance, it may be in- 
teresting to treat this problem in an unrestricted Hartree-Fock 
approximation, which can be used to systematically improve 
the variational wavefunctions used in this work. 

Edge states of surface alloys, such as the one discovered 
by Ast et al!^, seem to be promising candidates to realize 



the anisotropic Wigner crystal phase, since they combine ex- 
tremely strong Rashba SOC, a tunable Fermi energy, and 
screening due to the nearby metal. In a real system, however, 
disorder will inevitably play a major role. As a result, both 
the positional and orientational order are expected to be short- 
ranged. To detect the broken symmetry state on the surface, 
one can either resort to local probes (such as scanning tun- 
neling spectroscopy), or find a way to align the orientational 
domains, e.g. by an in-plane magnetic field. 
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Appendix A: Rashba particle in a box 

The arguments presented in this paper rely on the solu- 
tion of the quantum mechanical problem of a single parti- 
cle with Rashba SOC in a rectangular box of size L x x L y 
with infinite potential walls. While the corresponding prob- 
lem without spin-orbit coupling is trivial, with spin-orbit cou- 
pling the problem of boundary-condition matching with the 
multi-component wavefunction is highly non-trivial. The rea- 
son the problem with SOC is more difficult is that in this case, 
the Hamiltonian is no longer separable (i.e., it cannot be writ- 
ten as a sum of two commuting terms, one of which depends 
only on the x coordinate and the other on y). A circular well 
can be solved exactly 33 , thanks to its rotational invariance. 

In this Appendix, we combine several approaches to de- 
duce the asymptotic form of the ground state energy in the 
anisotropic, low-density limit, Eq. Q of the main text. We 
first solve the problem exactly in the L x — > oo limit (keeping 
L y fixed), and then provide an argument which yields the form 
of the leading corrections for finite L x . Finally, we present 
numerical results supporting the analytical arguments. 

1. Solution in the L x — > oo limit 

In the limit L x —> oo, the problem becomes translationally 
invariant in the x direction. The eigenfunctions then take the 
form 

ip(x,y) = e ik °> x tp(y), (Al) 

where ip, <p are two-component spinors. We choose coor- 
dinates such that the walls are at y — ±L y /2. Then <p{y) 
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satisfies the boundary conditions 



where 



<P(±Ly/2) = 



(A2) 



Seeking a solution with energy E — e, we get that the 



wavevector modulus k = ^/k% + k y satisfies 



(fc - k y 
2m 



(A3) 



Fixing k x , we find four allowed values of k y , which we denote 

by ±k y> ±: 



(A6) 



and a m , m (rji , 772 = ±) are coefficients which are determined 
by the boundary conditions. We may reduce the number of 
coefficients by using symmetry. Under reflection, y — >• —y 
the spinor ip(y) transforms as ip(y) —> a y ip(—y). Requiring 
that the wavefunctions are either even or odd under reflection 
gives 



k ±V2me) -k*. 



(A4) 



(A7) 



Note that k Vt ± can be imaginary. 

The wavefunction for the transverse motion takes the form 



<p(y) = a mV2 

r)!,r) 2 =± 



. , / p -iVl6 V2 /2 



(A5) 



Imposing the boundary condition, Eq. ( A2 1, on the wavefunc- 



tion in Eq. ( A5 I, for the even sector gives the following (im- 
plicit) equation for e: 



La 



ky-Ly 



6L 

T 



ky,-Ly 



(A8) 



r 



where k y ±, 6± are given by Eqs. ( A4i and ( |A6| . For the odd 
sector, we get an identical equation with cos replaced by sin. 
Equation (A8l determines the dispersion e(k x ). The disper- 
sion of the lowest subband, near k x — ko, is shown in Fig. [4] 
for widths k$L y /2 ranging from 10 to 20. 

The form of the low energy dispersion can be deduced from 
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Figure 4: Dispersion of the lowest subband for infinite strips of vary- 
ing width L y . The inset shows the dispersion minimum A vs. L y on 
a log-log plot, showing good agreement with A ~ 1/L* , 



Eq. ( A8 1 as follows. To lowest order in e and Sk x = k x — k , 
this equation depends on e, 5k x , and L y through the factors 

ky t ±L y : 



^y,± J -'y 



J (k ± %/2me^ - (k + 8k x ) 2 L y 



±2koV2me — 2koSk x L y . 
Therefore Eq. (|A8|> has the functional form 



(A9) 



F [2koV2meL 2 y , 2k 6k x L 2 y ) = 0, (A10) 

where F is some function of two variables. Close to the band 
minimum, the dispersion has the following form: 

2k V2meL 2 y = A 1 +A 2 [2k Sk x L 2 y - A 3 ] 2 , (All) 
where Ai 2.3 are dimensionless constants. Therefore 



e(Sk x ) 



1 



8mfcg 



^j+A 2 \ •2l,; )l )k, 



^\ L 2 



1 A\ A X A 2 



8mfc2 L y 



(6k x - sky 



(A12) 



where 8k* = A 3 /2koL y . We see the dispersion minimum A 

-4 
V 



scales as while the effective mass in the x direction is 



L independent. We confirm the relation A ~ 1/L y by solving 
Eq. ( A8 1 numerically, see inset of Fig. [4] 
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Figure 5: Numerical results for the ground state energy vs. aspect 
ratio L x /L y of a particle with Rashba spin orbit coupling in a rect- 
angular infinite potential well of size L x x L y , for several values of 
the density n (measured in units of 1/fco). We have set 2m = 1. 
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Figure 6: (a) Minimal ground state energy, e* , and (b) optimal aspect 
ratio, rj*, as a function of the density, on a log-log scale. The data 



is consistent with the analytical predictions: e 



,4/3 



-1/3 



and 77 



at low densities. 



2. Extension to Finite L x 



Next, we consider that L x is finite, still assuming that 
L x L y . In this case, the solution is complicated by mul- 
tiple reflections from the boundaries. Nevertheless, we can 
deduce the form of the ground state energy as a function of 
L x and L y as follows. For large L x , we assume that the solu- 
tion is largely composed of traveling wave states near the bot- 



tom of the dispersion given by Eq. ( A12 1, which are reflected 
back and forth from the two edges. Consider a right-moving 
wave with momentum k x = ko + 5k x . This state can only 
be reflected to the state k ' = k — 5k x + 25k*, where 8k* 



is defined below Eq. ( A12 1. Note that time reversal symmetry 
prohibits scattering to the other left-moving solution with mo- 
mentum —k$ — 5k x , since this state is the Kramer's partner of 
the original incoming wave. The eigenstates of the system are 
determined by the requirement that the phase acquired over 



one period is a multiple of 2tt: 

2 (6k x - 5k*) L x + ct> = 27tj, 



(A13) 



where j is an integer, and (f> = fa + fa is the sum of the 
two (unknown) phase shifts fa, fa associated with reflections 
from the two ends. The total phase shift is a function 5k x , 
L y and ko. However, in the limit 5k x L y — > 0, l/koL y —> 
0, we assume that we can replace cf> (5k x L y , l/koL y ) by a 
constant 4> (5k x L y — > 0, l/k L y — > 0) = fa. We then find 
that the ground state is given by 5k x — 5k* — (27rj m i n — 
fa)/2L x , where j m j n is an integer chosen to minimize the 
energy below. Inserting this into Eq. ( |A12| i, we finally get 

lT T s 1 A\ , A X A 2 (7rj min - fa/2) 2 



8mkl L y 



iLl 



(A14) 



Substituting n = \/(L x L y ) and 77 = L x /L y , we obtain 



Eq. (Ill 



3. Numerical solution 



In order to verify the assumptions that lead to Eq. ( A 14 1, we 
have numerically calculated the ground state wavefunction of 
a Rashba particle in a box. This is done using a generaliza- 
tion of the "plane wave decomposition" technique described 
in Refs. 1341351 The solution is written as a superposition of 
eigenstates of the free Hamiltonian, which are plane waves 
with wavevectors satisfying Eq. ( |A3[ ) for some value of e. The 
coefficients of the different plane waves are determined by re- 
quiring that both components of the wavefunction vanish on 
a set of points evenly distributed along the boundary (in this 
case, an L x x L y rectangle), and that the first component of 
the wavefunction spinor is equal to 1 at an arbitrary point in 
the interior. The sum of the squares of the wavefunction (the 
"tension") at a different set of points on the boundary is then 
calculated. The eigenvalues e are identified as the minima of 
the tension. The eigenvalues can be determined with an accu- 
racy of 1% or better. We have tested the technique by calcu- 
lating the eigenenergies of a Rashba particle in a circular well, 
and found excellent agreement with the exact results®! 

Fig. [5] shows the ground state energy vs. the aspect ratio 
V = L x /L y for a few values of the density 1/n = L x L y . 
The ground state energy is minimal for 77 ^ 1. We found this 
behavior for n < 0.08/cq; f° r higher densities (smaller box 
area), the minimum energy occurs at 77 = 1. In Fig. [6] we 
show the minimal ground state energy e* and the correspond- 
ing aspect ratio 77* for densities ranging from n/fcg = 8 x 10 -2 
to 4 x 10~ 3 . For low densities, the optimal ground state en- 
ergy and aspect ratio follow e* ~ n 4 / 3 and rf ~ n -1 / 3 , in 



agreement with Eqs. ( 10 1,( 111, and ( A14 1 



Appendix B: Hartree-Fock description of the ferromagnetic 
nematic 

Within the Hartree-Fock approximation, we replace the full 
Hamiltonian H, Eq.([T|i, by the following mean-field Hamilto- 
nian H: 
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2m 



2m 




(B5) 



where we kept only the leading order term in ^ . Next, we 

calculate the potential energy. Assuming that V (r) ~ r~ a 
at long distances, its Fourier transform has the following form 
for small momentum transfer: 



Figure 7: Constant e nerg y contours of the dispersion of the mean- 



field Hamiltonian, Eq. B2 around the minimum at k = k(,y. Numer- 



ical labels indicate energy values above the minimum, measured in 
units of fcg / (2m) . The dashed line indicates the circle k^ + ky — k% . 
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(Bl) 

where h is a spontaneous Zeeman field, to be determined 
self-consistently, and [i is the chemical potential. We pro- 
ceed by minimizing the expectation value of the full Hamil- 
tonian within the ground state of the mean-field Hamiltonian, 
Eq.pl). 

Let us focus on the case of an in-plane Zeeman field h. (We 
will later argue that an out-of-plane h is not energetically fa- 
vorable.) Without loss of generality, we assume that h = hx. 
Then, the lower-branch single particle dispersion obtained by 
diagonalizing Eq.(|Bl[) is 



£fc 



2m 




(B2) 



V(q)^V (Po-f3iq a - 2 ) 



(B6) 



for 2 < a < 4, where we have neglected higher order terms in 
q. Here /3 and j3i are constants, and Vq is defined in Eq.(|2]i. 
For a > 4, the leading order term goes as q 2 (as can be seen, 
e.g., from the fact that for a > 4 the second moment of the 
potential exists). We now compute the potential energy 



u = Je< f ^-^)> 



E 



n 3 



4+9,<r c M c k'- g ,«r' c *V') i ( B? ) 



k,k' ,g,<r,er' 

where Vl is the volume of the system, and we have introduced 
Ck, a , the annihilation operator of an electron with momentum 
k and spin a. The calculation can be simplified significantly in 
the limit of small mh/k" 2 ., in which the Fermi surface becomes 
parametrically eccentric. In that limit, only the dependence of 
V q on the momentum parallel to the major (long) axis of the 
Fermi surface is important. The result can be written as U — 
+ U^, where Uff and U^i are the interaction energies 
between same spins and opposite spins, respectively, given (to 
leading order in mh/k^) by 



The minimum of the dispersion is obtained for k = k y. In 
the low-density limit, only states close to the minimum are oc- 
cupied. We therefore expand the dispersion around the mini- 
mum, using polar coordinates, k x — k sin 8 and k y = k cos 9 
where k = k + 5k, to leading order in 6k and 6. This gives 



7.2 
ft 

~2m 



1 



2m 



1 



k 2 h 



4 (k 2 + fh) 



e 2 . (B3) 



The Fermi surface is thus approximately an ellipse centered 
around fcoy (see Fig|7]). Denoting the Fermi energy measured 

k 2 l 

relative to the dispersion minimum e = 5m + + M as 
ef, we get that the density and Fermi energy are related by 



2tt 



1 



2fcg 

mh 



(B4) 



The expectation value of the kinetic energy per particle in the 
ground state of H (measured relative to the band minimum, 

eo) is 



U n = A{p 1 V )n 
,(l3oVo)n 2 



U n = B- 



/■ 2 




a-1 



(B8) 



Here, A and B are dimensionless constants. We can now min- 
imize the total energy per particle, e = K + Uff + U^i, with 
respect to the variational parameter h. In the n — > limit, we 
find that h*, the optimal Zeeman field, is 



h* 



^(m^Vo^n 2 -^, 
P V n, 



a < 4 
a > 4. 



(B9) 



Inserting h* back into the expression for the total energy, 
we get that the ground state energy is 



iimVo^n 2 ^-^ 



a < 4 
a > 4. 



(BIO) 
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Finally, we argue that an out-of-plane Zeeman field, h z 7^ 
0, is not energetically favorable. First, note that for an in-plane 
field, in the low-density limit (n — > 0) the kinetic energy per 
particle is unaffected because the Zeeman field becomes par- 
allel to the spin-orbit field for all occupied states, allowing 
each particle to gain the Zeeman energy without changing its 
wavefunction. In contrast, in order for an electron to align 
its spin with an out-of-plane Zeeman field, its wavefunction 
must include hybridization with the excited band (with energy 
~ k%jm above the lower band). In the hybrized state, the up- 
per band is occupied with a probability ~ (m/i z /2/cq) 2 , lead- 



ing to a kinetic energy cost (fcg/2m) x (mh z /2k'^) 2 . For con- 
tact interactions, the potential energy per particle is given by 
E p = ^(U /n)(n 2 — ml), where m z s» nh z /(k /m 2 ) (here 
we use that for small h z , the spin of each electron obtains a 
z -component given by the ratio of h z to the spin-orbit field). 
Therefore the potential energy gain from z-polarization is ap- 
proximately ^Uon(mh z /kl) 2 . Compared with the kinetic en- 
ergy cost, we see that the potential energy gain includes an 
extra factor of n, indicating that the cost of polarizing in the 
z -direction overwhelms the benefits in the low-density limit. 
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